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Using a sample of 1.06 x 10® i// mesons collected with the BESIII detector at the BEPCII e+e" 
collider and Xcj mesons produced via radiative transitions from the 4>' , we report the first observation 
for ip' pK'^T," + c.c. (charge-conjugate), as well as improved measurements for the XcJ hyperon 
decays XcJ pK'^A + c.c. The branching fractions are measured to be B{tp' — > pK'^T,'^ + c.c) = 
(1.7±0.1±0.1)x 10-^ B{xcO -^pK+k + c.c.) = (13.2± 0.3± 1.0) x lO"*, B{xci -> p/f+A + c.c.) = 
(4.5 ± 0.2 ± 0.4) X 10"* and B{xc2 pK+A + c.c) = (8.4 ± 0.3 ± 0.6) x 10"*, where the first 
error is statistical, and the second is systematic. In the decay of Xco — > pK'^A + c.c, an anomalous 
enhancement near threshold is observed in the invariant mass distribution of pA + cc, which cannot 
be explained by phase space, and may be related to other such structures observed before. 
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PACS numbers: 13.25.Gv, 14.20.Jn, 14.40. Rt 

I. INTRODUCTION 

The study of hadronic decays of the cc states J/i/;, 
-0', and Xcj could provide valuable information on per- 
turbative QCD (pQCD) in the charmonium-mass regime 
and on the structure of charmonia. The color-octet 
mechanism (COM), which successfully described several 
decay patterns of the P-wave XcJ states [1], may be 
applicable to other Xcj decays. Measurements of Xcj 
hadronic decays may provide new input into COM and 
further assist in understanding the mechanisms of Xcj de- 
cays. Hadronic decays of charmonia below the DD mass 
threshold are also a good place to search for previously 
unknown meson states [2]. The BES Collaboration has 
previously reported observations of near-threshold struc- 
tures in baryon-antibaryon invariant-mass distributions 
in the radiative decay J — >■ ^pp [3] and the purely 
hadronic decay J/i/; — >■ pAK" ^ [4]. It has been sug- 
gested theoretically that these states may be observa- 
tions of baryonium [5] , or caused by final state interac- 
tions [6]. Studying the same decay modes in other char- 
monia may provide complementary information to im- 
prove the knowledge on these unexpected enhancements. 
It is also interesting to search for potential structures 
formed by A A and pT, pairs, which could assist in ex- 
tending the theoretical models. 

BESIII has gathered a large sample of e+e" — )■ ip' 
events, which leads to abundant production of Xcj states 
through radiative decays. This enables us to search for 
and study the hadronic decays of the Xcj states with high 
statistics. 



II. DETECTOR 

BEPCII [7] is a double-ring e+e^ collider that has a 
peak luminosity reaching about 6 x 10'^^ cm^^s"^ at a 
center of mass energy of 3770 MeV. The BESIII [7] detec- 
tor has a geometrical acceptance of 93% of Att and has 
four main components: (1) A small-cell, helium-based 
(40% He, 60% CgHg) main drift chamber (MDC) with 



t Throughout the text, inclusion of charge conjugate modes is im- 
phed if not stated otherwise. 



43 layers providing an average single-hit resolution of 135 
/j,m, and charged-particle momentum resolution in a 1 T 
magnetic field of 0.5% at 1 GeV/c. (2) An electromag- 
netic calorimeter (EMC) consisting of 6240 CsI(Tl) crys- 
tals in the cylindrical structure barrel and two endcaps. 
The energy resolution at 1.0 GeV is 2.5% (5%) in the 
barrel (endcaps), while the position resolution is 6 mm 
(9 mm) in the barrel (endcaps). (3) Particle Identifica- 
tion (PID) is provided by a time-of-flight system (TOF) 
constructed of 5-cm-thick plastic scintillators, with 176 
detectors of 2.4 m length in two layers in the barrel and 
96 fan-shaped detectors in the endcaps. The barrel (end- 
cap) time resolution of 80 ps (110 ps) provides 2<t K/tt 
separation for momenta up to ^ 1.0 GeV/c. (4) The 
muon system (MUC) consists of 1000 m^ of Resistive 
Plate Chambers (RPCs) in nine barrel and eight endcap 
layers and provides 2 cm position resolution. 

III. MONTE-CARLO SIMULATION 

Monte-Carlo (MC) simulation of the full detector is 
used to determine the detection efficiency of physics pro- 
cesses, optimize event selection criteria, and estimate 
backgrounds. The BESIII simulation program [8] pro- 
vides an event generator, contains the detector geometry 
description, and simulates the detector response and sig- 
nal digitization. Charmonium resonances, such as J/tp 
and ip' , are generated by KKMC [9, 10], which accounts 
for the effects of initial-state radiation and beam energy 
spread. The subsequent charmonium meson decays are 
produced with BesEvtCen [11, 12]. The detector geom- 
etry and material description and the transportation of 
the decay particles through the detector including inter- 
actions are handled by Geant4 [13]. 

IV. DATA ANALYSIS 
A. Event selection 

Candidate V' — > pK^^P and i/"' ^ IXcj ^ jpK+A 
events, with S*^ — > 7A and A — )• pn^ , are reconstructed 
using the following selection criteria. 

Charged tracks must have their point of closest ap- 
proach to the beamline within ±30 cm of the interac- 
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tion point in the beam direction {\Vz\ < 30 cm) and 
within 15 cm of the beamhne in the plane perpendic- 
ular to the beam {Vr < 15 cm), and must have the polar 
angle satisfying | cos 6*1 < 0.93. The time-of-flight and en- 
ergy loss dE/dx measurements are combined to calculate 
PID probabilities for pion, kaon, and proton/antiproton 
hypotheses, and each track is assigned a particle type cor- 
responding to the hypothesis with the highest confidence 
level (C.L.). For this analysis, four tracks identified as p, 
p, if"*", and 7r~ are required. To suppress backgrounds 
from fake tracks, the p and are constrained to the 
same vertex by vertex fitting, and are required to satisfy 
1141 < 10 cm and Vr <1 cm in the case oi^pK^ K modes, 
and the same procedure is applied for the respective an- 
tiparticle combinations in the charge-conjugate mode. 

Photon candidates are selected in the EMC by requir- 
ing a minimum energy deposition of 25 MeV within the 
barrel region | cos 6*1 < 0.8, and 50 MeV within the endcap 
regions of 0.86 < | cos 6*1 < 0.92. EMC cluster timing re- 
quirements suppress electronic noise and energy deposits 
unrelated to the event. 

A kinematic fit that enforces momentum and energy 
conservation (4C) is applied with the hypothesis i}}' — >■ 
"fppK^n" , where the p and tt^ are constrained by A 
decay vertex fitting. For the events with more than one 
photon candidate, the combination with the smallest xlc 
is retained for further analysis. 

A candidates are selected by requiring the invariant 
mass of p7r~ to be within 7 MeV/c^ of the mass of the A 
as given by the PDG [14], and this distribution is shown 
in Figure 1. E*^ candidates are formed by calculating the 
invariant mass of 7 and A candidates, and this is shown 
in Figure 2(a). 

After vetoing tp' pK^Yf' events by removing events 
where the 7 and A have an invariant mass within 15 
MeV/c^ of the S" mass [14], Xcj{J — 0, 1,2) signals are 
seen distinctively in the spectrum of recoil mass against 
the 7, as shown in Figure 2(b). 

B. Background studies 

For the measurements of Xcj ~^ pK^A, a sample of 
1.06 X 10^ inclusive ip' MC events are used to investigate 
possible backgrounds. The surviving events can be classi- 
fied mainly into three decay processes: (1) ip' — pK~^A, 
where a fake 7 is produced; (2) -0' — Tr^pX+A where 



one 7 from the 7r° decay escapes detection; and (3) the 
direct decay ip' — ^ 'ypK'^A having the same final topol- 
ogy with the signal, but not going through an intermedi- 
ate Xcj state. Accordingly, 2 x 10^ MC events for each 
of the three background processes are produced for fur- 
ther detailed studies. The same selection criteria are ap- 
plied to the exclusive MC samples, and the surviving 
events are normalized to 1.06 x lO'* total ip' MC events. 
For the normalization procedure, the branching fraction 
B = (1.00 ±0.14) x lO"'^ for 0' ^ pK+A is quoted in the 
PDG and the other two background modes have branch- 
ing fractions in the order of 10~^, which we roughly deter- 
mine from our actual data sample. Figure 3(a) presents 
the distributions of the recoil mass against the 7 for 
events that survive all cuts for the data and also for these 
background exclusive MC samples. 

A similar study is also done for the measurement of 
ip' —7- pX+E" using the three background modes above 
together with tp' — >■ jXcj ~^ jpK^A — > 'yppK'^ir^ , as 
shown in Figure 3(b). 

In addition, a 42.9 pb~^ data sample collected at 
3.65 GeV is used to investigate possible continuum back- 
grounds. Only 7 events survived inside the mass region of 
Xcj for the measurements of Xcj pK'^ A, and are found 
to be negligible. For ip' pi^T+E", 110 events from the 
continuum contribution must be subtracted after proper 
normalization according to the luminosities. 

C. Determination of branching fractions 

1. Number of ip' ^ pK'^TP events 

The decay mode ip' — >■ pK'^YP is observed for the first 
time, with the main background processes ip' — 7p_R'''"A, 
Ip' 7r"pi4r+A, Ip' pK+A and ip' -> -fXcj jpK+A. 
According to the studies in the previous section, the 
background shape can be described by a linear function, 
as shown in Figure 3(b). 

A maximum likelihood fit is applied to the spectrum 
of the invariant mass of the selected 7 and A, and we 
find a yield of 276 ± 21 events for the E" signal. The 
shape of the E° is obtained from MC simulation where 
the mass and width are fixed to the PDG values. The 
derived curves are shown in Figure 4, where dots with er- 
ror bars represent the data with continuum contribution 
subtracted. 
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FIG. 1. (Color online) The invariant-mass distributions of pn 
the A peak. 
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FIG. 2. Distributions of (a) the invariant masses of 7A and (b) the recoil mass against the 7 in decays of ^p' after vetoing 
^' pK+TP events. 



The detection efficiency for this process is determined 
to be 24.4% from MC simulation with a phase space 
model. The invariant mass spectra of pYP and YPK'^ 
are shown in Figure 5. 

2. Number of tp' — > 7Xcj — >■ jpK'^A events 

For the XcJ pK^A decays, obvious inconsistencies 
exist in the distributions of pK'^ and AK^ invariant mass 
between the phase space MC and data, as shown in Fig- 
ure 6, so the detection efficiencies for the decay modes 
V'' ~^ 7Xc0.ci,c2 jpK^A are determined by taking into 



account the dynamics of the decay. 

For each Xcj state, the allowed regions of M{pK^) 
versus M{AK~^) are divided into 25 x 25 areas of equal 
length (40 MeV/c^ for Xco and 48 MeV/c^ for xa and 
Xc2), and each area is tagged with an index ij. For each 
area the number of events N^^^^ for data and detection 
efficiency are determined individually. Then, the total 

number of events (A^cor) is calculated as iVcor — 
Samples of 5.5 x 10^ MC events are used to determine 
the detection efficiencies of each area for XcOi Xci, Xc2, 
respectively. 

The data belonging to XcO 1 Xci , and Xc2 are separated 
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FIG. 3. Comparison of data with exclusive MC samples for distributions of (a) the recoil mass against the 7 for — > jXcJ 
"fpK'^A and (b) the 7A invariant mass for ^' pK+Tp. The MC samples have been normalized to the total number of tp' 
events. In figure (a), the background from t]/ — >■ pK'^A events is too small to be visible. 




FIG. 4. (Color online) The shape of the E° signal as derived from MC simulations which had the mass and width fixed to the 
PDG values. The fit result is shown by the solid line with a linear background indicated by the dashed line. The data points 
with error bars show the data, where the continuum contribution has already been subtracted. 



using mass windows on the distribution of recoil mass 
against the detected 7 of 3.35-3.48, 3.49-3.53, and 3.53- 
3.59 GeV/c^, respectively. When extracting -/V^ata' 
background has been subtracted using exclusive MC sam- 
ples according to the results of background studies. The 
calculated total numbers of events A^cor are listed in Ta- 
ble I. 



TABLE I. The total numbers of events Ncor for each Xcj ^ 
pK A are derived from A'^cor ~ "Sij—f^^. A^crror is the prop- 
agated error. 



Modes 


A"cor 


error 


XcO 


8642.7 


201.3 


Xcl 


2824.0 


112.6 


Xc2 


4961.0 


154.4 
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3. Calculation of branching fractions 

The branching fraction of -0' — is calculated 
with 

g _ ^obs 

where N^i is the total number of ijj' events, which is mea- 
sured to be 106 X 10® with an uncertainty of 0.81% [15]; 
the branching fractions (63.9±0.5)% for Ba^ptt and 100% 
for B^a^^/^ are taken from the PDG [14]; iVobs means the 
observed number of signals derived from the fit and e is 
the detection efficiency from MC simulation. 

The branching fractions for each XcQ,c1.c2 — ^ pK~^A are 
calculated similarly with 



where the branching fractions of the Xcj states ((9.68 ± 
0.31)%, (9.2±0.4)% and (8.72±0.34)% for 6(V'' ^ 7Xco), 
B{'4)' — )■ 7Xci) and B{tp' — > 7Xc2), respectively) are taken 
from the PDG [14]. 



D. Near-threshold structure 

The large discrepancies between the data and phase 
space MC samples in Figure 6 imply that intermediate 



states exist in the decays of x^j pK^A. Possible struc- 
tures are observed in the Dalitz plots shown in Figure 7, 
and particularly for the Xco , it seems that there is a struc- 
ture in the near-threshold region of M{pA) reflected by 
the anomalous enhancement in the top right corner of 
the Dalitz plot. 

Figure 8(a) shows the invariant-mass distribution of 
pA for XcO ~^ pK~^A, where the dashed line denotes the 
phase space distribution that has been normalized to the 
signal yield and the dots present efficiencies in each bin. 
Evident discrepancies are seen near the threshold region. 
Due to insufficient statistics, in this analysis a simple fit 
with a Breit-Wigner function to this region is done with- 
out considering quantum mechanical interference. The 
fit curve for the near-threshold structure is depicted in 
Figure 8(b), where the distribution of M{pA) has been 
corrected by the detector efficiency. The structure can 
be fit well with a weighted Breit-Wigner function of the 
form 

^^^^ (AP-MSr-M^T^ 

where q is the anti-proton momentum in the pA rest 
frame, k is the kaon momentum in the XcQ rest frame, 
L [L') denotes the orbital angular momentum between 
the antiproton and A (between the kaon and pA). On 
the basis of conservation on J^, in the decays of XcOj 
"L + L' = even number" can be inferred, and therefore 
the only possible spin-parity combinations are — 0~ , 
1+, 2^, • • • . Because the structure is near the pA thresh- 
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FIG. 6. Invariant mass spectra of pif"*" and AK~^ for (a, b) Xco, (c, d) Xci and (e, f) Xc2- The dots are the data, and the 
hatched regions show the distribution of MC events generated according to a phase space model. Potential intermediate states, 
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FIG. 7. Dalitz plots of M^{pK^) versus A4^{AK^) for (a) XcO, (b) Xci and (c) Xc2- A concentration of events in the upper 
right corner shows an enhancement at the pA threshold. 



old, the relative orbital angular momentum between the 
antiproton and A is most likely 0. Therefore, = 0~ is 
used in the fitting process which gives M = 2.053 ±0.013 
GeV /c^ and T = 292 ± 14 MeV for the Brcit-Wigner 
mass and width parameters. A shape of the phase space 
MC is added to describe the background in the fitting, 
which is shown as the dashed line in Figure 8(b). 



For tp' — >■ pK^'S'-', the invariant-mass spectrum of 
M(j?E°) was shown in Figure 5(a). In this channel, there 
may be similar structures close to the threshold, but 
there is a large uncertainty due to the relatively small 
sample size. 



V. SYSTEMATIC UNCERTAINTIES 

The main contributions to the systematic uncertainties 
in the measurements of the branching fractions originate 
primarily from the tracking, PID, photon reconstruction, 
kinematic fit, branching fractions of intermediate states, 
total number of ip' events, and the fitting procedure. The 
results are summarized in Table II. 

The tracking efficiency for MC simulated events is 
found to agree with the data within 1% for each charged 
track coming from a primary vertex from analyses of 
J/tp — > K*K and J/ip — pp7r+7r~ events. For each track 
from A (or A), the uncertainty is also 1% according to a 
study of very clean J/V' — >■ pK~^A events. 

The candidates for the selected final states require 
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FIG. 8. (Color online) (a) Invariant-mass distribution of pA for XcO — >■ pK^A, where the dashed line denotes the phase space 
distribution that has been normalized to the signal yield. The histogram shows the data and dots present the efficiency curve, 
(b) Fit result to a Breit-Wigner function with — 0^ after acceptance correction. The dashed line describes the background 
shape from phase space MC events. 



tracks to be identified as p, p, or tt^. Compar- 
ing data and MC event samples for J/t/j — > p/f+A and 
J/ip —i' K*K, the difference between MC and data for 
the particle identification efficiency was found to be 2% 
for the antiproton, 1% for the proton and kaon, and neg- 
ligible for charged pions. 

The difference in the reconstruction efficiency between 
the data and MC is about 1% per photon [16]. 

To estimate the uncertainty from kinematic fitting, the 
kinematic fitting efficiency is studied using events of "0' 
7XcO ~^ ^ppTT^TT^ and the difference between data and 
MC is found to be 2.8%. 

Uncertainties due to the mass window requirement for 
the A signal are studied with the control sample tp' — >■ 
pK^A. The efficiency difference between data and MC 
is obtained to be 0.4%. 

Uncertainties in the fitting procedure are obtained 
by varying fit intervals and changing the linear back- 
ground shape to a 2"^^ order Chebyshev polynomial. It 
contributes a 1.8% uncertainty to the measurement of 

The uncertainty on the total number of events was 
found to be 0.81% by studying inclusive hadronic ■0' de- 
cays [15]. 

Uncertainties due to the branching fractions of ■0' — >■ 
7XcJ are 3.2%, 4.3% and 3.9% for each XcO, Xa and Xc2, 



respectively [14]. The uncertainty due to the branching 
fraction of A ^ p7r" is 0.8% [14]. 

Uncertainties due to the numbers of areas in the pro- 
cedure of calculating total numbers of events for t/j' — > 
7Xc,/ •ypK'^A are shown as "2D Binning" in Table II. 
Detection efficiencies are assumed to be constant within 
each of these 25 x 25 sub-areas (see section IV C 2), and 
as a check, we varied the number of areas. Besides the 
original 25 x 25 binning, three other divisions (20 x 20, 
30 X 30, 35 X 35) were tried, and the largest differences 
among them are taken into account as the systematic 
uncertainty due to the binning. 

Uncertainties from the mass window requirements of 
XcO J Xci and Xc2, obtained by changing the Xcj selection 
window, are shown as item "Mass Window" in Table II, 
and are small compared to other errors. 

A possible A polarization in the decays of Xcj might 
affect detection efficiencies and yield different results. 
With our limited statistics, it was not possible to mea- 
sure the polarization of the A in fine bins of the Dalitz 
plot for each Xcj state, but an overall measurement of 
the A polarization P was done for each Xcj state that 
yielded P 0.04 ± 0.07 for XcO, -0.17 ± 0.12 for Xci, 
and 0.22 ± 0.09 for Xc2- Subsequently, new samples of 
MC events were then generated with the A having this 
polarization P, so that the decay distributions are given 
by 1 -I- aPcosO, where Q is the angle between the A 
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flight direction in the XcJ rest frame and the tt direction 
in the A rest frame, and a is the weak decay parameter 
for the A. The difference in efflciencies with respect to 
that of phase space MC samples are taken as a systematic 
uncertainty. 

The total systematic uncertainty is obtained by sum- 
ming up uncertainties contributed from all individual 
sources in quadrature. 



VI. RESULTS AND DISCUSSION 

We observe the decay mode ij]' — > pX+E" + c.c. for the 
first time and improve the measurements for the decays 
of Xc.j pK^A + c.c, using 106 x 10^ ?/'' events col- 
lected with BESIII detector at the BEPCII colhder. The 
branching fractions are listed in Table III. 

For the pK~^A + c.c. final state in the decays of Xco, 
an anomalous enhancement is observed in the invariant- 
mass distribution of pA + c.c, which could correspond to 
the structure observed in the decay J/ip — pAK~ [4]. 
It is of great interest that the structure is located very 
close to the mass threshold of pA + c.c., and this may be 
accounted for as a quasibound dibaryon state or as an 
enhancement due to a final-state interaction, or simply 
as an interference effect of high-mass TV* and A*. Our 
new measurements may aid in the theory of charmonia 
decays, and also be a guide in the calculation of decay 
modes into strangeness dibaryon systems. A more precise 



measurement is expected with larger statistics in future 
BESIII running. 
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TABLE II. Systematic uncertainties in percent (%) 
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Xcj pK+A 
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Tracking 
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PID 
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Photon Recon. 
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0.81 


2D Binning 




1.3 
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